THERMODYNAMICS AND UNIVERSALITY FOR MEAN FIELD 

QUANTUM SPIN GLASSES 



NICHOLAS CRAWFORD 
Department of Mathematics, University of California at Los Angeles 

ABSTRACT. We study aspects of the thermodynamics of quantum versions of spin glasses. By 
means of the Lie-Trotter formula for exponential sums of operators, we adapt methods used to 
analyze classical spin glass models to answer analogous questions about quantum models. 



1. Introduction 

Classical spin glass models have seen a flurry of activity over the last few years, see [2, 4, 9, 
10, 11, 13, 14, 15, 20, 21, 22] to name a few; in particular [2, 9, 10, 22] all consider aspects of 
the (generalized) Sherrington-Kirkpatrick model, a mean-field model in which the interactions 
between spins are mediated by an independent collection of Gaussian random variables. In con- 
trast, though physicists have considered both short and long range quantum spin glass models 
for quite sometime, see [3, 8, 18, 19], there are few rigorous mathematical results; we mention 
[5], which provides a proof that the quenched free energy of certain short ranged quantum spin 
glasses exists. 

Here we extend classical spin glass results to quantum models in two directions. First, using 
the ideas of [11], we demonstrate the existence of the quenched free energy of the Sherrington- 
Kirkpatrick spin glass with a transverse external field. Next, under conditions made precise below, 
we give a complementary result which shows that a large class of quantum spin glasses, including 
the transverse S-K model, satisfies universality. By this we mean that as the size of the system 
goes to infinity, the asymptotics of the free energy of the system do not depend on the type of 
disorder used to define model. This latter result is based on the work of [4]. 

One may view this paper as an attempt to adapt the methods of classical spin systems to the 
analysis of various quantum models. Guerra's interpolation scheme [9, 10, 11] and the Gaussian 
integration-by-parts formula are ubiquitous tools in the classical setting. A major theme of the 
present work is that through the systematic use of the Lie-Trotter product formula (i.e. that 
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for any pair of operators A, B) one may extend the interpolation scheme and integration-by-parts 
formula to quantum systems in useful ways. 

For concreteness, the remainder of the introduction and the following section use the language 
of the spin- 1 /2 representation of su(2) to describe quantum spin systems, though a number of 
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our methods apply in larger generality. In this setting, one describes each particle using a two 
dimensional Hilbert space C 2 along with a representation of su(2) generated by the triplet of 
Pauli operators S = (S (x \S^,S^). 

To represent an N -particle system, we introduce the tensor product Yn = ®f=i C 2 , one 
factor for each particle, along with a sequence (Sj)^^ of N copies of the Pauli vector S, where 
Sj acts on the j'th factor. The particles interact by means of the Gibbs-Boltzmann operator 
e ~/3H N assoc i a ted to the Hamiltonian Hn- Here, Hn is a self adjoint operator acting on Vn, 
typically a polynomial in the iV-tuple of spin operators (Sj)jL. v For example, the Hamiltonian 
of the simplest non-trivial quantum spin system, the transverse Ising model, is described by 



i,j=l j=l 

where A > 0. 

Once we specify the Hamiltonian, statistical quantities of the system may be defined. The 
partition function and free energy of a quantum spin system are defined via the trace of the 
Gibbs-Boltzmann operator as 

Z N (J3) =Tr(V^) 

Self adjoint operators on Vat replace functions as observables of the system and the thermal 
average of an observable A is defined as 

Tr(Ae _/3Wjv ) 

With this formalism we present a few examples of spin glasses of particular interest. Tradi- 
tionally the modeling of any spin glass necessitates the introduction of disordered interactions 
between spin. In general, and will be the case here, the interactions are i.i.d. The basic example, 
the transverse S-K model, has a Hamiltonian defined by 

N N 



2^fN .4^ ' J 1 3 3 

1,1 = 1 7 = 1 



*J = 1 3= 

A more complicated class of models, the quantum Heisenberg spin glasses, are described by one 
of the Hamiltonians 

1 N 
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where a, 7 € 
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There is a sharp and interesting contrast in the scope of application of our results. Unlike the 
quantum Heisenberg spin glasses, applying the Lie-Trotter expansion to the transverse S-K model 
has the added benefit of allowing a natural path integral representation, known as the Feynman- 
Kac representation, of statistical quantities like the free energy. More precisely, the Feynman- 
Kac representation allows the expression of the partition function as a measure on cadlag paths 
with state space the classical Ising spin configurations on N sites, {—1, +1}^. For this and other 
technical reasons we can only resolve the existence of the free energy for the one quantum model. 

The remainder of the paper is structured as follows. Section 2 introduces notation and states 
the main theorems of the paper. Section 3 presents quantum generalizations of [4] which allow 
us to prove universality and control the fluctuations of the free energy for a class of quantum 
spin glasses which include the spin- 1 ^ quantum Heisenberg spin glasses. In addition, we adapt 
the techniques of [11] to prove exponential decay of the fluctuations of models with Gaussian 
disorder. Finally, Section 4 details the existence of the pressure of the transverse S-K model. 



2. Results 

We begin by giving a bit of notation to be used below. We consider a collection of random 
variables, where / is in some finite indexed set 3 and denote E [•] and P (•) integration with respect 
to this collection. Examples of index sets that we have in mind include the collection of all r- 
tuples of sites in a system of N particles. In the simplest case the index set 3 consists of all pairs 
(i, j) (or more generally some subset of pairs). Unless otherwise specified, the variables £/ are 
assumed to be i.i.d. according to some fixed random variable £ satisfying the conditions 

E[£] = 0, E[f] =1, E[|e| 3 ] <oo. 

When explicitly considering Gaussian environments we denote the random variables by gj. 
Let S_ represent the iV-tuple of Pauli vectors and consider the Hamiltonian 

where each Xj(S_) is a self-adjoint polynomial in the spin operators, i.e. Xj(S_) = Xj(S_). We 
define the associated partition function and quenched 'pressure' by 

Z N (J3, = Tr (e^) ; a N (J3, £) = E [log Z N (J3, £)] ■ 

Note that for convenience we have omitted the minus sign from the expression for the Gibbs- 



Boltzman operator. For any operator A we denote its operator norm by \\A\\ := sup^ GV 
With this notation we have the following result: 



(Av,Av) 
(v,v) 



Theorem 2.1 Let £ be random variable with mean 0, variance 1, and E |£| 3 < oo and let g 

be a standard normal random variable. Let 3 n index the interactions between particles at the 
system size N. Suppose 

' ||X 7 GS)|| 3 = o(A0. 
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Then for any (3 € M, 



as N tends to infinity. Moreover 

1 



E 



v logZ(/3,e)--a^v(/3,e) 



< 



o(l) 

AT3 



Remark 2.2 In the quantum Heisenberg Hamiltonians, the norm of each summand is bounded 
by -jL and the number of pairs is of the order N 2 . As a result, 



£ ||x 7 gs)|| 3 = o{m\ 



ie3 N 

and the theorem is immediately applicable. 

The error bounds get better if one assumes the random variable £ has higher moments. At the 
extreme end, we consider the case of fluctuations for Gaussian environments: 

Proposition 2.3 Let g be a standard normal random variable. Then 



F(\log Z N ((3,g)-a N ((3,g)\ > u) < 2 e V 

Remark 2.4 Note the generality with which these results are stated. In the introduction we ad- 
vertised their application to mean field models. However, choosing the operators Xj and index 
sets 3 jv appropriately allows application to a wide variety of systems. 

Next we take up the more subtle question of the existence of the free energy for mean field 
models. As mentioned above, we specialize to the transverse S-K model: 

N N 



i,j=l j=l 

In light of Theorem 2.1, we have assumed the interactions are Gaussian 
Theorem 2.5 Let (3, A > be fixed. Then 

-1 



lim E 

jv^oo (5N 



logTr (e-PH N (*y 



exists and is finite. Moreover, there exists a K > so that the following concentration property 
holds: 



1 



TV 



logTr (e"^" (A) ) + — E 



N 



log Tr ( e 



-0Hn{\) 



> u) < 2e 



—N- 



Remark 2.6 For readability, we assume here that the disordered portion of the Hamiltonian in 
(??) is a two body interaction. However analogous arguments allow one to treat p-spin models 
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where we replace the disordered portion by 

oo r 

-««•<«<•>) -"E^E*,.* IK' 

r=l i\,...i r k=l 

Here {g^...^} is a collection of independent standard Gaussian random variables and Y^T a rQ r 
is even, convex and continuous as a function of q on [— 1, 1]. See [2] and [11]. 



3. Universality and Fluctuations 

Before proving Theorem 2.1, we adapt the methods of [4] so as to apply them to a wide variety 
of quantum spin systems. It turns out that the line of argument given there is robust enough to 
be followed in the quantum case, though the calculations must be adjusted to accommodate the 
non-commutative setting. 

In the general setup, we consider a collection of self adjoint operators {Xi}f =l and Ho defined 
on some finite dimensional Hilbert space. Let 

d 



Z(/?,0=Tr(e^«)) 
a(/U)=E[logZ(/?,£)] 



denote the Hamiltonian, partition function, and quenched 'pressure' of a system. We define the 
thermal average, Duhammel two point function, and the three point function for the operators 
A, B, and C as follows: 

v _Tr(Ae^) 

/l 1 Tr (Ae^^Be^^O) du 

(AB) - WJ) ~ 

(A, B, C) _ 1 . 

For the convenience of the reader, we recall a generalized integration-by-parts lemma proved 
in [4]. The left hand side of (3.1) is well known to be zero if the randomness is Gaussian. This 
fact will play a role in what follows. 

Lemma 3.1 Let £ be a real valued random variable such that E [|£| 3 ] < 00 an d E [£] = 0. Let 
F : M — > R be twice continuously differentiate with ||-F"||oo = sup^gjj < oo. Then 



|E[£F(0] -E [^ 2 ] E [F'(0]\ < ^"llooE [^r 
We first provide an expansion of the quenched pressure: 



(3.1) 
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Lemma 3.2 



d(3 

where \0{(3 2 )\ < (3 2 . 



(3E 



.1=1 



+ 90(/3 2 )E [|^| 3 ] (^IIX, 



Remark 3.3 There is no correction in this formula when £ is Gaussian since the integration-by- 
parts formula is exact in this case . 

Proof. To begin the derivation, we have 



i=l 



Let us define Hi(z) = J2j^i ZjXj + zXi + ^Hq, denote the corresponding Gibbs state by (-)^ 

(z) 

and define the function Fi(z) = (X,^ . Then 

To expand Equation (3), we calculate the first and second derivatives of Fj. Applying the Lie- 
Trotter expansion, i.e. that e A+B = lim fc e A/fc e S/fe , with ^ = and B = - /^JQ, 
we find that the first and second derivatives of Fj take the form: 



F'M =(3 
F>i(z) =(3 2 



{X^p-^X^) 
2 (X i} X l} X^ - 3 (X t , X0 W (Xi) W + 2 (W*))' 



In order to prove the lemma we must bound F". To this end we claim that for any a\ , . . . , a n > 
so that ai H + a n = 1 and any self-adjoint operators X, 7i 

|Tr (le aiW • • • Xe a " H ) \ < \\X\\ n Tr (e H ) (3.2) 

k ■ 

Indeed, we may assume by continuity that a,j = for some m > and some sequence of 
positive integers (kj) summing to 2 m . 

2k 



For any sequence of operators (BM^, 




< 



2k 

n Tr ([W) 



k\ 2k 



(3.3) 



j'=i 



This can be seen as a specific case of the general method of chessboard estimates, Theorem 4.1 
of [7]. Applying (3.3) to the left hand side of (3.2) with Bj G {Xe^, e^} we have 



|Tr (Xe aiH ■ ■ ■ Xe anH ) \ < Tr 



Xer^X 



Tr (e n ) 



, 2" l -n 

ri\ 2 m 



(3.4) 
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Another application of (3.3) implies 



Tr 



Using this bound on the right hand side of (3.4) and letting m pass to infinity proves the bound 
(3.2). 

It follows that ||-F/'||oc < 6/3 2 ||Xi|| 3 . Recalling that E [£ 2 ] = 1, Lemma 3.1 in conjunction 
with the previous calculations imply 



d(3 



(3E 



1=1 



<9/3 2 E[|£| 3 ] 



vi=l 



□ 



Next we use this expansion to compare the quenched 'pressure' for £ to that of a Gaussian 
environment g where g is a standard normal. 

Lemma 3.4 Let and {gi\ be collections of i.i.d random variables distributed according to 
£ and g respectively. For any (3 G M, 



\a((3,0-a((3,g)\ < 9|/3| 3 E [|£| 3 ] £ 



\i=l / 

Proof. The similarity between the proof of this lemma and that of Proposition 7 [4] means that we 
will be extremely brief. Consider the interpolating partition function and corresponding quenched 
'pressure' defined by 

Z( s ,t- s) =e v / 5(Eti6X l )+v / ^(Eti^ l )+w„ 
a {t) (s) =Elog Z(s,t- s) 

respectively. We have 

<*(*>(*) = a(Vt,£); a«(0) = a(A<?)- 
Lemma 3.2 along with independence between the two environments implies that for all s £ [0, t] 
we have 



ds 



<9^E[|e| 3 ] l^WXi 



,i=l 



For (5 > 0, if we let t = j3 2 and integrate this inequality the result follows. For [3 < we instead 



consider the environments — £, — g. 



□ 



Next we attend to the fluctuations of the 'pressure' determined by the random environment 

(6): 

Lemma 3.5 77zere exists some universal constant c > so that 

E [| logZ(/3,£) - a(/3,£)| 3 ] < cE [|£| 3 ] /3 3 Vd ^ ||^|| 3 ) . 
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Proof. Consider the filtration = a{£i . . . k > 1 determined by the sequence of indepen- 
dent random variables Let 



A, :=E[logZ(P,0\fi] -E[\ogZ((3,0\Fi-i] 



We have 



logZ(/3,£)-a(/U) = E A - 
Burkholder's martingale inequality implies the existence of a universal constant d so that 



E 



i=i 



< c 'e(^a 



vi=l 



To bound the increment A,, consider the partition function 

Zi(j3,0=Ti(eP**) 
where ipi = 4>i(0 '■= — Since Zi((3, £) is independent of £j, 



Ai =E 



log 



— E 



log 



Zi(j3,0 



i-l 



We use this identity to estimate Aj. 
We claim that 

Indeed, the Lie-Trotter formula implies 

(ok 
T e 2* e 2fe 



< e /?lfcHI*ill. 



lim Tr TTe 2k+1 e^'e 2 fc+1 
fc^oo \ 1 1 



Since 



e 2fe 



< e 2*= 



Inequality (??) applied to the righthand side for k finite gives 



u=i 



from which our claim follows. 
From this we estimate: 



I Ai| </?||X;||(|&|+E|&|). 
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Therefore, 

/ d \ § / d 

E|logZ(/3,0 -a(/3,0| 3 < CE E A H ^ E M 2]E (I6I + E l&l) 5 



vi=l 



vi=l 



<c/? 3 E|&| 3 ^ E W X * 



□ 



Proof of Theorem 2.1. Our first theorem now follows as an application of the above machinery: 
the first statement is an application of Lemma 3.4, while the second follows from Lemma 3.5. □ 

Finally, we consider fluctuations in Gaussian environments: 

Proof of Proposition 2.3. Let Zp(t) be defined as the auxiliary partition function given by two 
independent collections of Gaussian disorder and g^ 2 \ 



*(*)=Tr(e 



with t an interpolation parameter varying between and 1. Let Ej denote the average with respect 
the random variables for j = 1,2. Given any s£R, let 



Y(t) = exp( S E 2 logZ /3 (t)); <f>(t) = logEi [Y(t)] . 

We note that 

<p N (l) - <Pn(0) = logE[exp( S (logZ(/?,<7) -a{P,g)))] . 
In order to estimate this difference, consider 

d 

Y(t)E 2 ' 



(3.5) 



2Ei [T(t)] 



Ei 



where the notation (-) t represents the Gibbs state induced by the interpolating Hamiltonian. A 
calculation involving the Lie-Trotter expansion and the Gaussian version of the integration-by- 
parts formula implies 



f' {t) - wmf 1 



i/t 



i=l 



so that 



„2o2W WY-W 2 

W(t)\ < s(j hi^M , 

Using Equation (3.6) and the fact that e^l < e x + e~ x we have 



(3.6) 



exp(\s\\logZ((3,g)-a(f3,g)\)<2e 
Finally, applying Markov's inequality we have 



>(\\ogZ(J3,g)-a(p,g)\ >u)<2e 



^ 2 /3 2 £f = l ll^ill 2 
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for any seR. Optimizing over s concludes the lemma. □ 



4. Existence of the Pressure 

Recalling the su(2) formalism, we choose a preferred basis for Vjv consisting of tensor prod- 
ucts of eigenvectors for the operators {S^}. Denoting the eigenvector for which corre- 
sponds to the eigenvalue +1 by |+) and the eigenvector for which corresponds to the eigen- 
value — 1 by |— ), we may identify this preferred basis with classical Ising spin configurations 
a G {—1, +1}^. For each a, we denote the corresponding basis vector by \a). 

The proof of Theorem 2.5 proceeds in two steps. The first step consists of a concentration 
estimate following essentially the same argument as that of Proposition 2.3. Widening the scope 
beyond the transverse Ising spin glass, for this first step we consider quantum Hamiltonians of 
the form 

n N :=n% s {s^) + np (4.i) 

We assume that only 7i% s involves Gaussian disorder and that the deterministic operator Tiff 
takes a sufficiently nice form so as to admit a Feynman-Kac representation in terms of the basis 
of eigenvectors for the operators. By this we mean that 

(a\exp(-un d ^)\a) > 

for all u > and all spin configurations a, a. Assuming that all off diagonal matrix elements of 
Ti^ are negative gives a necessary and sufficient condition which guarantees the existence of a 
Feynman-Kac representation. As mentioned in Section 2, for our treatment Tiff takes the form 

N 



i,j=l 

where the collection {gij} is assumed to be i.i.d. standard normal, though more general interac- 

(z) 

tions involving the {Sj } may be considered. 

To illustrate the significance of the Feynman-Kac representation, recall that by adding a suit- 
ably chosen diagonal matrix to H d ^ we can force the rows of the matrix representation of Tiff 
to sum to 0, which implies that the (modified) Hamiltonian Hf} 1 is the generator of a continuous 
time Markov chain with state space {— 1, +1}^. 

Consider any matrix element (o"|e~ MWjv \a). Expanding the exponential for finite k using the 
Lie-Trotter formula with A = —uTiff(S^) and B = —uTijf and inserting the complete or- 

A B 

thogonal set {\cr')} between each factor e^e^ and then passing to the limit in k, it is not difficult 
to see that 

(a\e- uH »\a)= [ Jo -^H)^ dl/(<7 ) 

where v is the induced Markov chain measure and 

H% s (a) = (a\H% s (S^)\a) 
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Moreover, we can view the augmentation by the diagonal matrix as introducing a weighting to 
each spin configuration which corresponds to the amount of time the process spends in each state 
along its trajectory. Thus, the original Hamiltonian yields an un-normalized measure on paths 
taking values in { — 1, +1}^. Let us also note, as is the case in the transverse S-K model, that the 
diagonal matrix is constant. 

Returning to our original Hamiltonian (4.1), the upshot is that via the Feynman-Kac transfor- 
mation we may represent the partition function Zn associated to the Hamiltonian TLn by 

Z N =Tr(e-^) = I e J " ^N Z ^^ Mu)aAu)du du(a} 

where v is a measure on spin configuration paths determined by 7i A ^ and U is the space of 
cadlag paths taking values in [—1,1]. We use this representation implicitly throughout Section 
4. In the case of the transverse Ising model, one may check that the induced measure is in fact 
(proportional to) the Markov chain measure defined by starting from an initial configuration and 
evolving in time via a spin flip process determined by flipping the spin value at each site according 
to the arrivals of independent Poisson processes of rate A. In order to prove convergence of the 
quenched pressure we shall need to consider restricted partition functions Z N defined via the 
Feynman-Kac transformation by 

_Lio g Z# = llog J J" ^^^ 9ijai{u)aj{u)du du(a) 

where A = An is some deterministic (i.e. not depending on Gaussian disorder) subset of f2. 
For any pair of spin configurations a and a, let 



1 N 



a denote the spin path {a(u) : u G [0, (3]} and 

h8 1 V 



V 1 

Cn(z) = / -7= Yl 9ij^i(u)^j(u)du 
Jo 2VN i j=l 



denote the Gaussian random variable defined by the classical Gibbs weight of the spin path a. An 
easy calculation shows that 

N f 13 f@ 

E [OvGiKat®] = - \ \ R 2 (a(u),a(s))duds 
4 Jo Jo 

where the function R is the classical overlap for Ising spin configurations. 

Lemma 4.1 Let (3 > Obe given. Fix any Feynman-Kac representable deterministic Hamiltonian 
H%'. Then 



llog^-lE [\ogZ N ] 



>u) <2exp( - 



2u 2 
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Proof. With the appropriate modifications, the method of proof of Proposition 2.3 may be fol- 
lowed here: Let Z N (t) be defined as the auxiliary partition function given by two independent 
collections of Gaussian disorder g^> and g( 2 \ 

Z N {t) = [ e ^N^ + ^N^du(a) 

J A 

with t an interpolation parameter varying between and 1 and £W the Gaussian corresponding 
to g®. Let Y N (t) and <f> N (t) be defined in terms of Z N (t) as in (??) and (??). 

Replacing the corresponding quantities appearing in the proof of Proposition 2.3 we have, 
using the Feyman-Kac representation, 



s 2 



M) ~ 2 El [Y N (t)} 

Ei Y N (t) / du{a)du(d) - / / i? 2 

[ JA N xA N V 4 Jo JO / 

where a) is the truncated 'Gibbs weight' corresponding to the event An. Thus 

The bound now follows as in the proof of Proposition 2.3. 



□ 



Remark 4.2 More generally, the p-spin models may also be treated via the method employed 
here, though the bound stated in the lemma must be modified slightly. 

Unfortunately the use of the Feynman-Kac transformation alone does not allow our method 
to go through. In particular, we were unable to treat the quantum system with deterministic 
quadratic couplings in the x and y directions: the ferromagnetic version does permit a Feynman- 
Kac representation but the interaction is convex, which turns out to have exactly the wrong sign 
in the expression for the derivative of the interpolating 'pressure'. 

The only natural example that we found amenable to our method is the transverse field Ising 
model. Notice that the deterministic portion of this particular Hamiltonian is linear, which simpli- 
fies the interpolation scheme that we employ to analyze the thermodynamics at different system 
sizes. For inverse temperature f3 and transverse field strength A > 0, we refer to the partition 
function of this model by Zn((3, A) and denote pn{(3, A) = — jj log Zn(/3, A). 

Lemma 4.3 Let j3, A > be fixed. Then 

lim E\p N (P,X)]=p(P,\) 

N~*co 

exists. 

Proof. This proof, like that of the previous lemma, relies on the proof of an analogous statement 
in [1 1]. The main idea is to partition the space of paths into subsets on which we may control the 
time correlated self overlap R(a(u),a(s)). 
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To this end, let e, 5 > be fixed, where for convenience we assume j € N. For any function 



g : [0, (3] x [0, (3] -> [-1, 1] we define the event A g (S, e) by 



A,(5,e) ={a : g(iS,jS) < R(a(i5),a(j5)) < g(i5,j5)+eVi,j < | - 1, 

<2eV( U ,s) G [0,/3) x [0,/3)}. 

Let S'a(e) be the set of functions which are constant on [jS, (j + 1)5) x [kS, (k + 1)<5) for j, k G 
{0, ... | - 1} and take values in {ie : i e i] n N}. We define the event 

A = A(S,e) = U geSs (e)A g (5,e). 

Observe that though A definitely does not cover the full sample space Q, it is enough to prove 
convergence of the truncated pressure 



p£(AA) = -llog^e^W(£). 



More precisely suppose e = e^,5 = 5^. We claim that if —e^logS^ is sufficiently large as 

N -> oo 

lim E[p N (P,\)-p£(P,\)] =0. 

Jy — »oo 

Let AJ\fis(<r) denote the total number of jumps made by the spin path a in the time interval 
[iS, (i+l)S]. To determine how large to take — e log S with N, let A* = {a: m.&x i< p_ 1 AA/^a) > 

^}. Then ^ c C A* so that 

J^exp (Cjv(g))di/(g) X4 it exp(Cjy(oi))di/(a) 
/4exp(Civ(2))di/(£) - J A cexp(CAr(<i))dz/(o;) 
By Jensen's inequality, 

J Am exp(( N (a))dis(a) 



E 



log 1 + 



< 



f A c exp (Cjv(2)) dz/ (^) 

log | 1 + / E 

By the Cauchy-Schwarz inequality, 

exp(0v(V)) 



exp(CAr(q:))d^(a) 
Lc exp(Cjv(oi))d^(o:) 



di/(o) . (4.2) 



E 



/ As exp(6v(a))di/(£} 

< E [exp (2(n(o))]^ E 

<E[exp (2C N (a))]^E 



(j A cexp(( N (a))dv(afj 
eXP (^)L CiV(£)dl/(£) 



c\-l 
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1 Acdu 



The last line follows from Jensen's inequality applied with respect to the path measure -^J^y 
Since 2(jv(V) is a Gaussian random variable with variance N R 2 (a(u), a(s))duds, 

N ^ rP 



E[exp(2Cjv(£))] =exp 
Similarly, after a short calculation we have 
( 2 



o jo 



R 2 (a(u),a(s))duds 



E 



exp 



( N (a) dv(a) 




As a result of the preceding estimates 

exp(0v(V))<M£) 



f E 

J A* 



J AC ex.-p(( N (a))dv(a) 



R 2 (a(u),a(s))duds 



^ v{A if ) 
du(a) < e 2 



dv{q)dv{a) 

^M 2 



Standard calculations imply that for all e small enough and 5 < e 2 , 



< I p CiVelog«5 

S 6 



(4.3) 



for some universal constant C > 0. Requiring —CelogS > ft 2 H — -^log<5, putting estimates 
together and taking the appropriate limits proves our claim. Note that these conditions can be 

arranged, for example, by letting ejy = N~ 1/4 and 6n = e~ N /2 and taking iV large. 

Thus we are reduced to showing that the mean of the truncated pressure E [p^{P, A)] con- 
verges. For convenience of exposition we first consider subsequences of the form 7V fe = N n k for 
some Nq, n G N. For any k, we may view a G &,N k as an n-tuple of spin paths (a^, . . . , £ ( - n * > ) 
so that cjW G ^iv fc _i- hi order to compare thermodynamics at consecutive system sizes, let us 
define the interpolating Hamiltonian 

n 

c Nk (t,a) = vtc Nk (a) + VT^Y,dlM { ^ 

1=1 

where (jy involve disorder couplings which are mutually independent and independent 

from the couplings in (a) . In addition, we introduce the partition function 

Z&(t) = [ exp(( N (t,a)) 
J A 

Let 

A g = A Nk (5,e,g)={a = (a^,...,a^) : a® € A^S^g)} 
Obviously A g C A g . Therefore 



-logZ£(t)<-logZ£(t) 



(4.4) 
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For any g G Ss(e) consider 



After bit of work employing the integration-by-parts formula we arrive at 



+ -E 



where (-)f 3 corresponds to the truncated Gibbs weight determined by the Hamiltonian £jv fc (t, a) 

and (-)f corresponds to the product Gibbs weight determined by an independent pair of spin 
paths a, a_ and Hamiltonian (jy k (t,a) + 0v fe We stress that corresponding terms in this pair 
Hamiltonian involve the same realizations of disorder. From the definition of A g the first term 
can be bounded by (3 2 e in absolute value. As / is convex, the latter term is less than or equal to 
zero. Thus, by evaluating (f>$ at zero and one and applying (4.4) we have 



P 1 P f p 

/ R 2 (a{u),a(s))duds - - V / / R 2 (a ( - l) (u), a(s))dud 
Jo n i Jo Jo 

/ / R 2 (a(u),a(s))duds--Y, [ [ i?V (u),* (0 (*)) ) 
Jo Jo n i Jo Jo / 



— E 



< 



N, 



-E 



k-l 



log 2^ +/3 2 6 



for any g G S^e). 

Next, by Lemma 4.1, we have 



^ log (/?, A) - i-E [log (0, A) 
Setting u = e, there exists a set S and a universal constant D > so that 



> u) < 2exp ( "^iVfc ) • 



P(«S) > 1 - exp 



2e^ 

' f3 2 n 



and so that on S 



-1 



— log Z% ((3, A) < — log Zg^ 09, A) + D(3 2 e 



N, 



k-i 



for all g G 55(e). 

Now let us choose (5 = exp(— 2 (p e )2 ) an ^ e 
Borel-Cantelli lemma implies that there exists D > so that 



N k 4 . Then a standard application of the 



-1 



-1 



— log Z^ k (J3, A) > — log (J3, A) + D/3 2 AT- 1/4 for some g G 5 5(JVfc) (iV- 1/4 ) i.o. ) = 0. 
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Therefore, for large enough (with the choice of e and 8 depending on as above) we 
have 



E 



,A{8,e,g) 



(P, A) 



> <r D ^ lS ' 

ges s (e) 



; / 4 



> e 



i(^A) 



where we have used the inequality xf > k n (Xa Xi 



for n > 1 and Xj > 0. 



By (4) we have 



-E 



log (/?, A) - log Z N ^((3, A) 



Therefore the truncated sequence of 'pressures' — log Z^ Nk ' tNk \(3, A) is nearly decreasing 
a.s. It follows that the truncated thermodynamic limits exist almost surely. 

Let Xk = — jj- log Z N ^ k (/?, A). Since {E [X^]} is uniformly bounded, the limit is finite a.s. 
Another application of the concentration inequality (4.1) along with the Borel-Cantelli lemma 
implies that this limit is a non-random constant f((3,\). Further, it is now a small matter to 
prove convergence of the truncated quenched averages along these subsequences. The bound 
(4.1) implies that X^ are uniformly integrable, i.e. 







lim supE [|X fe |1| Xfe |> K ] = 0, 



N 



v 4 



so we may conclude E [X^] — > /(/3, A) as well. Finally, as — €N k log 5iy k = the arguments 
given above imply that lim^ E \p^ k (/3, A)] = f(/3, A). Another application of the concentration 
inequality then implies that PN k ((3, A) — > f(/3, A) a.s. 

It remains to prove uniqueness of this limit and then convergence along arbitrary subsequences. 
However these statements follow from fairly standard arguments based on the above results. □ 
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